A quantitative measurement of spatial order in 
ventricular fibrillation 



P.V. Bayly E.E. Johnson P.D. Wolf 
H.S. Greenside W.M. Smith 
R.E. Ideker 

From the Engineering Research Center 
for Emerging Cardiovascular Technology 
School of Engineering 
Duke University 
Durham, North Carolina 27708 
and the Departments of Pathology and Medicine 
Duke University Medical Center 
Durham, North Carolina 27710 

Supported in part by the National Institutes of Health 
research grants HL-28429, HL-44066, and HL-33637 

and the National Science Foundation 
Engineering Research Center Crant CDR-8622201 



February 9, 2008 



Abstract 



Introduction: The degree of organization in ventricular fibrillation (VF) is 
not known. As an objective measurement of spatial order, spatial correla- 
tion functions and their characteristic lengths were estimated from epicardial 
electrograms of pigs in VF. 

Methods: VF was induced by premature stimulation in 5 pigs. Electrograms 
were simultaneously recorded with a 22 x 23 array of unipolar electrodes 
spaced 1.12 mm apart. Data were obtained by sampling the signals at 2000 
Hz for 20 minutes immediately after the initiation of VF. Correlations be- 
tween all pairs of signals were computed at various times. Correlation lengths 
were estimated from the decay of average correlation as a function of elec- 
trode separation. 

Results: The correlation length of the VF in pigs was found to be approxi- 
mately 4-10 mm, varying as fibrillation progressed. The degree of correlation 
decreased in the first 4 seconds after fibrillation then increased over the next 
minute. 

Conclusions: The correlation length is much smaller than the scale of the 
heart, suggesting that many independent regions of activity exist on the 
epicardium at any one time. On the other hand, the correlation length is 4 to 
10 times the interelectrode spacing, indicating that some coherence is present. 
These results imply that the heart behaves during VF as a high-dimensional, 
but not random, system involving many spatial degrees of freedom, which 
may explain the lack of convergence of fractal dimension estimates reported 
in the literature. Changes in the correlation length also suggest that VF 
reorganizes slightly in the first minute after an initial breakdown in structure. 
Key words: spatial correlation, correlation length, organization, dynamics, 
ventricular fibrillation, sudden death. 
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1 Introduction 



The ventricular myocardium is a nonlinear, excitable system which can ap- 
pear almost totally synchronized in sinus rhythm, or completely disorganized 
in fibrillation. Fifty years ago Wiggers described four stages of ventricular 
fibrillation; his classification relies strongly on qualitative descriptions of the 
degree of spatial organization such as "undulatory contractions" (stage 1), 
"convulsive incoordination" (stage 2), "tremulous incoordination" (stage 3), 
and "feeble contraction wavelets" (stage 4) Reentrant spatial structures 
such as wavelets or rotors are thought to provide the underlying mechanism 
of VF. Rotating activation patterns have been observed in experiments in 
animals and in isolated myocardium p|, 0, they have also been widely 

studied in mathematical models of excitable media [P, 0, §]. In general, 
greater spatial order is associated with fewer, larger structures. In this paper 
we explore an objective measure of the amount of spatial organization: the 
spatial correlation length. 

Many investigators have examined the organization of activation in 
the heart by detecting potentials at a number of locations in the myocardium 



with an array of spatially distributed electrodes ^ |T^. Organization 
and spatial structure have been documented by following the progress of 
individual waves of activation. Isochronal maps, which depict the position 
of the wave front as time progresses, or isopotential maps, which show lines 
of constant voltage on the surface of the heart at a single instant, illustrate 
the size and type of waveforms which occur during VF. Isochronal maps are 
more commonly used than isopotential plots, since one map can show the 
position of the wave front at many instants of time. However, objective, 
statistical descriptions of activity are also needed to compare with results 
from simulations and from other experiments. 

The question of spatial organization is related to the widely speculated 
relationship between nonlinear dynamics, chaos, and fibrillation [^1], 13 



The mathematical term "chaos" refers specifically to the temporally aperi- 
odic, unpredictable behavior of some nonlinear deterministic systems charac- 
terized by broad-band power spectra and extreme sensitivity to initial condi- 
tions. A necessary condition for chaos is at least one positive Lyapunov expo- 
nent, which describes the exponential separation of nearby orbits. The con- 
cepts of strange attractors and chaos have illuminated the behavior of simple, 
few degree-of-freedom, nonlinear systems. The methods of low-dimensional 
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chaos also apply to the nonlinear dynamics of multi-degree-of-freedom sys- 



tems which involve very few modes [Tj] such as an oscillating fluid system 

Temporally complicated, but spatially organized, behavior of extended 
systems may often be well modeled by low- dimensional chaos. In these cases 
one may obtain useful quantities, like Lyapunov exponents, and analyze sim- 
plified models which retain the basic mechanisms of instability and disor- 
der. Spatially incoherent dynamics in an extended system are inherently 
high-dimensional, meaning that many independent variables are required to 
describe its state. These cases are impractical to treat as chaotic |]T6[. Be- 
cause of this difference, determining the degree of spatial organization of VF 
is important in developing useful models of fibrillation and finding effective 
clinical countermeasures. 

In other fields, notably fluid mechanics, a numerical measure of the 
characteristic size of a spatial process has frequently been extracted: the cor- 
relation length. There is more than one definition of the correlation length; 
the quantity. A, used here describes the rate of decay of the spatial corre- 
lation function as a function of distance. It has been calculated for many 
experiments, and has been interpreted as the characteristic length scale of 
the typical spatial structure in the experiment: the diameter of the average 
wavelet, eddy, or rotor, for example |T8|, |I9|, ^ 



While VF has been the subject of much study and more speculation, 
an objective, quantitative description of its spatial organization has not been 
put forward. Clear experimental results are needed to compare with the 



burgeoning supply of data from numerical models |2^, In this study, the 



correlation length of ventricular fibrillation has been calculated for several 
intervals, at various times after the initiation of VF, in each of five pig hearts. 
The application of the method and the experimental procedure are explained 
in more detail below. 



2 Methods 



2.1 Spatial correlation function 

The spatial correlation function, R{dx.), is a scalar function describing the 
average correlation between all pairs of measurements separated by the vector 
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(ix. For a general, real- valued function of time and space, g{t,x.), the spatial 
correlation function is 1231: 



Ridx) = E{g{t,^)git,x + dx)}, (1) 

where E denotes an ensemble average. For a stationary, ergodic process, 
the value of the spatial correlation function is obtained from the following 
integral. Integration is over the time interval — T to T and over the spatial 
domain S. 

/2(c/x) = lim — / / (yfft, x) (yfft, X + (ix) (it (ix. (2) 

T-»oo 2T Js J-T 

The spatial correlation function is completely analogous to the commonly 
encountered temporal autocorrelation function, 

A{t) =\im ^ r f{t)f{t + T)dt. (3) 

Sampling the continuous process, g{t, x), with a 2-dimensional array of equally 
spaced sensors leads to a set of discrete measurements. If P time samples are 
taken from each sensor of an M x array, one obtains the values: 



Cijk = g{tk, Xi,), z = 1, 2, M, J = 1, 2, N,k = 1,2,..., P, (4) 

where tj. is the time of the k^^ sample and Xjj is the location of the sensor in 
the i'^^ row and j^'^ column of the sensor array. For a regularly spaced array 
with a separation of h units between adjacent sensors, Xjj = [h{i — l),h (j — 
1)]^, and dx = [hAi,hAjf. 

The spatial correlation function of this discrete data set is ideally 
approximated by the sum: 

MNP 

R{Ai, Aj) = ^^^^ ^ ^1 ^ CijkC^,+Ai)ij+Aj)k- (5) 

Here Q{Ai, Aj) is the total number of simultaneous samples which are sepa- 
rated in space by the the vector [Ai, Aj]'^. The quantitities Ai and Aj may 
vary from — (M — 1) to {M — 1) and from —{N — 1) to (A^ — 1), respectively. 

The resulting estimates of spatial correlation are only statistically 
valid for values of Ai and Aj which allow an adequate number of correlations 
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to be averaged; for example, in an M x array there are only two sensors 
separated by a displacement of [(M— 1), (A^— 1)]^, but MN sensors separated 
by the displacement [0,0]"^. Accordingly, correlations will not be presented 
for separations larger than 3/4 of the shorter side of the rectangular array. 

Just as the temporal autocorrelation function of a harmonic signal is 
a sinusoidal curve with the same period, the spatial correlation function of a 
travelling plane wave is a sinusoidal surface with the same wavelength (Figure 
1). The spatial correlation function of uncorrelated random noise is close to 
zero everywhere except at the origin (Figure 2). In many real, physical 
systems (turbulent fluids, capillary waves), the envelope of the correlation 
function does not remain constant, as in the perfectly ordered plane wave, 
nor does it vanish as in the case of the random noise. Instead, it usually 
decays at some characteristic rate. For small Icixl the decay can be assumed 
to be exponential, as: 

i?(H) oc e-N/\ (6) 

To illustrate how we extract the correlation length. Figure 3a shows 
a typical decaying correlation, _R, plotted as a function of the vector sepa- 
ration between sensors, [dx, dyY . The solid curve of Figure 3b is obtained 
by averaging the surface of Figure 3a over all angles to obtain a plot of R 
vs |(ix|. The decay of the correlation curve is described by the parameter. A, 
obtained from a least-squares exponential curve fit also shown in the Figure 
3b. Two numbers can be extracted from the original curve of correlation vs 
separation: Lg is the length needed for the correlation to decay to 1/e, and 
Lq is the value of |(ix| at which the correlation first drops to zero. may be 
used as a check on the value of A. 

Although the exponential curve is not always the optimal fit, it is 
a simple, straightforward, repeatable, and widely used approach. Values of 
A vary quite closely with values of Lg (see Figure 3 and Table 1). It is 
impractical to use Lg or Lq if not all curves decay to 1/e or 0. 

A wide variety of processes have been characterized by their correla- 
tion lengths. Hohenberg and Shraiman in reference used the correlation 
length. A, to describe the degree of coherence of one of the simplest nu- 
merical models of disorder in a spatially extended nonlinear system: the 
1-dimensional Kuramoto-Sivashinsky equation. Tufillaro and colleagues in 
T9| studied the correlation length of a system of small surface waves in wa- 



ter as a forcing parameter was varied. Coullet and coworkers in [PTI] , recently 
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used A as a quantitative measure of spatial organization in numerical studies 
of the complex Ginzburg-Landau equation, a 2-D system which exhibits spi- 
ral patterns. They show that for their model, the spatial correlation length 
corresponds specifically to the mean distance between spiral "defects" , which 
is presumably also the average diameter of the structures Although no 
systematic evidence is presented of such a specific meaning for the correlation 
length of fibrillation, we interpret it as roughly the average length of a patch 
of coherent activity. 



2.2 Experiment 

Fibrillation data were obtained from 5 anesthetized pigs, each weighing ap- 
proximately 25-27 kg. The animals were anesthetized with pentothal 5%, 
and respiration was supported by a mechanical ventilator. Anesthesia and 
fluid levels were maintained with an intravenous line. A line from the femoral 
artery allowed the continuous monitoring of hemodynamic status and blood 
gas values. The animals received doses of the drugs ketamine and acepro- 
mazine before anesthesia was started, and intravenous succinylcholine was 
used to induce and maintain skeletal muscle paralysis. Normal blood gas, 
electrolyte, and body temperature levels were maintained. 

The chest was opened using a median sternotomy, and the heart was 
exposed and suspended in a pericardial cradle. The electrode plaque was 
sutured to the epicardium of the right or left ventricle. The electrode plaque 
is a plastic block containing 524 Ag-AgCl electrodes, 506 of which form a 
22 X 23 array. The electrodes in the grid are spaced 1.12 mm apart, providing 
total dimensions of 2.35 x 2.46 cm. A photograph of the array is provided 
in Figure 4. The electrode spacing was chosen to avoid significant spatial 
aliasing while allowing the plaque to cover as large an area as possible. A 
previous study has provided evidence that very little spatial aliasing error is 
introduced with 1.12 mm spacing pB| . 

VF was induced by applying a premature stimulus to the ventricular 
muscle. In this procedure the muscle was paced with relatively small, regular, 
current stimuli (SI stimuh) apphed to the myocardium through a stainless 
steel electrode. The stimulus electrode was sutured to the epicardium about 
1 cm from the electrode array. The SI stimuli were spaced 350 ms apart and 
each pulse delivered current equal to twice the diastolic theshold. After 10 
equally spaced SI pulses were applied, a larger (50 mA) current stimulus, S2, 
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was applied through another stainless steel electrode located 1 to 2 cm away. 
The S2 stimulus was applied prematurely, while some of the tissue was still in 
an inexcitable state from the previous SI beat. The large, premature stimulus 
is believed to induce rotors as it intersects the partially refractory region p|. 
This method is somewhat analogous to the procedure used by Davidenko and 
colleagues to initiate rotating waves in isolated sheep myocardium 0, and 
to the initial conditions used by Barkley and others in numerical models of 
other excitable media [^. Reference ||26[ describes the experiment in more 
technical detail. 

Simultaneous time series of the electrical potential at each of the elec- 
trodes were recorded with 14-bit precision at a sampling rate of 2000 Hz, 
using a 528-channel mapping system [^. The signals were prefiltered with 



a 500 Hz low-pass filter to prevent temporal aliasing and with a 0.5 Hz high- 
pass filter to eliminate low- frequency artifact. These frequency cutoffs are 
well removed from the characteristic frequency of VF in pigs, which is usu- 



ally near 10 Hz |2^. Potentials were measured relative to a reference voltage 
from the animal's left leg. 

For each experiment, 20 minutes of data were recorded immediately 
after VF was induced. Subsets of these data were taken from the first few 
seconds of fibrillation, and after 15 seconds, 30 seconds, 1 minute, 2 minutes, 
5 minutes, and 10 minutes of VF. Several spatial correlation functions were 
calculated using one-second blocks of data taken from each time interval. 
The one-second block size was chosen to ensure stationarity in time while 
including several cycles of activation; in several test cases the time interval 
was varied from 0.1 to 1.0 seconds, with only slight quantitative effects. 



3 Results 

A typical time series recorded at a single electrode just after the initiation 
of fibrillation is shown in Figure 5. Figure 6 displays a sequence of spatial 
patterns recorded by the electrode grid in the first second after the initiation 
of VF. The spatial oscillations do not appear to be smaller than the size of 
the array. Two seconds later the structure seems to have broken up, as the 
sequence in Figure 7 indicates. By 60 seconds, however, larger structures 
have reappeared (Figure 8). 

Spatial correlation estimates track these changes in organization con- 
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sistently. Figure 9 shows the spatial correlation computed from one-second 
blocks of data in the first, second, third and sixty-first seconds after the initi- 
ation of VF in one heart (experiment #3). In the first second the correlation 
surface falls off relatively slowly from the peak at the origin (Figure 9a). As 
disorganization increases (more structures of smaller size) the fall-off becomes 
more rapid (Figures 9b, 9c) and the surfaces approach the delta function of 
Figure 2. However, at 60 seconds the fall-off is noticeably less sharp (Figure 
9d). 

To extract a correlation length from the surface of spatial correlation 
the correlation estimate was averaged over all directions. Curves of average 
correlation vs. magnitude of separation are shown in Figure 10; these curves 
were obtained by averaging the spatial correlation surfaces of Figure 9. 

The curves in Figure 10 decay roughly exponentially, and were fitted 
using a least squares algorithm. The constant. A, for each curve is tabulated 
in Table 1. The length needed for the correlation to decay to a value of 
1/e, Le, and the value of the first zero-crossing of the curve, Lq, are also 
presented. It is apparent that A, Lg and Lq vary together. We can also 
see that the qualitatively more organized behavior has a longer correlation 
length. 



Time 


A 




Lo 


1 sec 


11.4 


10.2 


> 15.0 


2 sec 


6.3 


6.0 


10.8 


3 sec 


5.0 


4.6 


9.5 


60 sec 


8.7 


7.9 


12.7 



Table 1 : Characteristic lengths (in mm) obtained from the correlation curves 
of Figure 9 (experiment #3): the correlation length or constant of exponential 
spatial decay. A; the length to decay to 1/e, L^; and the first zero-crossing, 
Lo. 

Figures 11 and 12 show spatial correlation vs. vector separation and 
spatial correlation vs. the magnitude of separation for an episode of fibril- 
lation in a different heart (experiment #4). Although the progress of the 
arrhythmia is different, the characteristic correlation curves of Figure 12 are 
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qualitatively similar to those in Figure 10. The differences in the correlation 
length capture the quantitative differences in spatial organization. 

Similar analyses were performed for each of the five experiments and 
correlation lengths were calculated for data at various times during VF. Table 
2 provides a summary of the correlation lengths calculated at several different 
times, for each experiment. The correlation lengths fall mostly between 4 mm 
and 10 mm. 



Time 


Exp #1 


Exp #2 


Exp #3 


Exp #4 


Exp #5 


1 sec 


11.3 


5.9 


11.4 


6.1 


3.8 


3 sec 


4.4 


3.3 


5.0 


5.2 


3.7 


15 sec 


4.8 


3.7 


6.0 


6.0 


4.1 


30 sec 


4.8 


4.8 


7.6 


4.7 


3.7 


60 sec 


5.0 


5.4 


8.7 


8.1 


5.0 


120 sec 


6.0 


5.3 


6.7 


8.4 


4.6 


300 sec 


5.4 


4.8 


5.3 


6.7 


4.2 


600 sec 


5.7 


4.9 


5.7 


6.4 


4.1 



Table 2: Correlation lengths, A, (in mm) calculated from VF data during the 
progression of fibrillation in 5 experiments. Each value was calculated from 
a one-second block of data recorded at the indicated time after the initiation 
of VF. 

In this study it was found that the correlation length of VF decreased 
over the first few seconds but often increased as fibrillation progressed; in each 
experiment the correlation length after 3 seconds was shorter than both the 
correlation length after 1 second and the correlation length after 60 seconds. 
These differences were statistically significant with p < 0.05. The mean 
correlation lengths (over the 5 experiments) as VF progressed are shown in 
Figure 13. 

4 Discussion 
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4.1 Interpretation of results 

The correlation length is a widely used measurement of spatial organization 
which can be used to estimate the degree of spatial organization of physical 
processes. Longer correlation lengths are associated with larger structures 
and apparently more organized behavior. This study suggests that corre- 
lation lengths may be extracted from experimental epicardial electrograms 
and that those lengths may be used to compare different stages and types of 
ventricular fibrillation. 

Moe and his co-workers postulated in 1964 that the size and num- 
ber of wavelets in the "wandering wavelet" model of ventricular fibrillation 
was critical to the development of an episode of fibrillation ^9\. The broad 



interpretation of spatial correlation is not restricted to a certain physical pro- 
cess, nor does it explicitly provide evidence for a particular mechanism (such 
as wandering wavelets, or spiral waves). The physical mechanism should, 
however, be consistent with the specific numerical values of the correlation 
length. 

The initial breakup of organization during fibrillation was first ob- 
served years ago. It is described by Wiggers as the transition between "un- 
dulatory" motion and "convulsive incoordination" [Q. The apparent reorga- 
nization of VF after 30-60 seconds is not well known, but is consistent with 
the greater low-frequency content of typical temporal power spectra observed 
as fibrillation progresses ^ . Such restructuring may be due to a coales- 



cence of small structures into larger ones as time passes, perhaps related to 
the changing parameters of the myocardium. Worley and colleagues noted 
decreased activation rates in epicardial electrograms during the first few min- 
utes of VF in dogs; changes in activation rate on the endocardium were much 



smaller These effects were attributed to the ischemic condition of the 
heart. Numerical models have shown that even for fixed parameters, coales- 
cence of disorganized behavior into structures such as stable rotors can occur 



Can observed correlation lengths illuminate any aspect of the relation- 
ship between low-dimensional chaos and fibrillation? The correlation length 
of VF in pigs was found to vary between 4 and 10 mm. Simple scaling argu- 
ments can be used to establish a rough measure of the global complexity of 
fibrillation. If the correlation length corresponds to the length of the average 
coherent structure on the epicardium, each such region would have an area 
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of 16-100 mm^. The surface area of the ventricular epicardium in a typical 
pig heart may be about 7500 mm^. The ratio of these two areas implies that 
at a given instant, there exist 75 to 500 approximately independent areas 
of activity on the epicardium. This rough estimate of the number of coher- 
ent regions, although probably an upper bound (see Section |4^ ), provides 
evidence that VF is too complicated to be modeled as a low-dimensional, 
few-degree-of-freedom, system. Such evidence could explain the lack of con- 
vergence of the correlation dimension which has been reported by previous 
authors [|^, |3^. Our results suggest that treating VF as low-dimensional 
chaos is unlikely to succeed in models or in clinical strategies. 

On the other hand these results strongly support the claim, recently 
made in a study by Damle and co-workers [^, that VF in dogs is not random. 
Damle et al. showed that the directions of activation recorded by adjacent 
electrodes (separated by 2.5 mm) were related. This is consistent with our 
observations because the adjacent electrodes were closer than a correlation 
length apart. Damle et al. also document the decline in the amount of 
spatial "linking" over the first 5 seconds of VF. They did not examine longer 
episodes of VF. The finite correlation lengths we observed also provide direct 
evidence that VF is not a spatially random process. Instead, both studies 
indicate that VF is a high-dimensional process, but with a measurable degree 
of coherence. 



4.2 Limitations 

The estimate of 75 to 500 independent regions calculated above might be 
interpreted as a generous upper bound on the number of activation fronts 
which exist simultaneously on the epicardium. The size of the electrode ar- 
ray, which covers perhaps 5-10% of the area of the ventricular epicardium, 
prevents the observation of high correlations between events further apart 
than the length of the array. Thus, for behavior with a characteristic wave- 
length larger than the size of the array, the estimate of correlation length 
is biased downwards. Also, our statistic measures correlations between si- 
multaneous events, so two sites exhibiting similar activity at different times 
will appear under-correlated. Finally, the choice of the unsealed correlation 
length A, rather than say 2A, to represent the length of a "correlated" region 
represents a conservative and somewhat arbitrary decision. 

Averaging over the area of the array and over one-second blocks of 
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time may diminish the contributions of local or short-lived organized pro- 
cesses. However, the spatial correlation function is independent of any tem- 
poral variability during the averaging time. For example, as long as their 
spatial properties are the same a stationary spiral wave will have the same 
correlation length as a drifting spiral, though their local temporal signals will 
be quite different. As another example, if an identical, temporally random, 
signal were recorded by each electrode, the electrodes would still be perfectly 
correlated. 

No evidence has been presented of a direct relationship between the 
correlation length and any other specific quantity measured during VF. Some 
of the conclusions drawn in this paper are based on analogies to other dis- 
tributed systems: notably fluid experiments and numerical simulations. Im- 
proved experiments and methods of analysis may allow the calculation of the 
mean distance between reentrant loops, or a similar quantity, to firmly estab- 
lish the absolute meaning of the correlation length in fibrillation. A possible 
length scale for comparison is the width of the region of electromotive activity 
during fibrillation, which was recently reported by Witkowski and colleagues 



34 1 to be between 0.5 and 1.0 mm in dogs; this quantity is roughly an order 



of magnitude smaller than our estimated correlation lengths. 

Although our measurements were taken from a 2-dimensional array, 
VF occurs in a 3-dimensional medium. It has been shown that well-defined 
3-D structures can produce apparently spatially complex patterns of activity 
on a 2-D surface ||35| . Therefore it is quite possible that even though the 2-D 
patterns may be well characterized, the underlying 3-D process is not. 



4.3 Implications 

It is potentially extremely useful to measure quantities which are associated 
with organization. If the break-up or coalescence of reentrant waves is part 
of the mechanism of fibrillation or defibrillation, then a measure of the co- 
herence of VF must be important to the clinician or the electrophysiologist. 
Drugs may be compared by studying whether, or how much, they increase 
or decrease spatial organization. Does more organization make it more or 
less difficult to defibrillate? Does the optimal defibrillation strategy depend 
to the spatial nature of the arrhythmia? Is a high degree of organization 
associated with a tendency to spontaneously defibrillate? Is the size of a 
rotor related to the size of the smallest heart that will sustain fibrillation? 
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It is possible to answer these questions only if it is possible to compare un- 
biased measures of spatial structure. The correlation length is a repeatable 
and objective measure of relative organization. The description of spatial 
organization in VF by such a quantity is a useful step toward understanding 
and effectively treating fibrillation. 
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Figures 



Figure 1. 

The spatial correlation function of artificial data which are periodic in time 
and 2-dimensional space: f{x,y,t) = sin{4t — 2x — y). a) The spatial pattern, 
f{x, y, to), observed at one instant in time (a snapshot of the traveling wave), 
b) The correlation, R, as a function of displacement, or spatial lag, [dx, dyY . 
The data consisted of 100 time samples at each "sensor" of a 15 x 15 spatial 
array. 

Figure 2. 

Same as Figure 1. except that data which was generated by a UNIX random 
number generator: f{x,y,t) = randomO. a) Random instantaneous spatial 
pattern, b) Correlation vs vector displacement: the spatial correlation is 
approximately a delta-function (zero everywhere except at the origin). 

Figure 3. 

a) A typical surface of correlation vs. separation vector (spatial lag), b) The 
corresponding curve of correlation vs. magnitude of sensor separation. Also 
shown are the characteristic constants A, Lg, and Lq. Data was from 1000 
samples of VF at each sensor of a 22 x 23 array. 

Figure 4. 

Photograph of the electrode array. Silver/silver chloride electrodes are ar- 
ranged in a 22 X 23 grid with 1.12 mm inter-electrode spacing. 

Figure 5. 

A representative time series of VF recorded at one electrode, one second 
after fibrillation was induced by a premature (S2) stimulus. The sampling 
frequency was 2000 samples/sec. 

Figure 6. 

Spatial patterns of electrode potentials recorded by the 22 x 23 array during 
the passage of an activation front 1 second after the beginning of VF. Darkest 
pixels represent potentials of less than -2 millivolts; lighter pixels represent 
progressively higher potentials up to +2 millivolts; potentials of +2 millivolts 
and higher correspond to the lightest grey level. ("Bad" electrodes are re- 
sponsible for some isolated anomalous pixels) . Panels represent "snapshots" 
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of the array at a) t=1.230 sec. b) t=1.250 sec. c) t=1.270 sec. d) t=1.290 
sec. 

Figure 7. 

The spatial patterns of electrode potentials recorded by the array 3 seconds 
after the beginning of VF. a) t=3.210 sec. b) t=3.230 sec. c) t=3.250 sec. 
d) t=3.270 sec. 

Figure 8. 

Patterns of electrode potentials 1 minute after VF was induced, a) t=61.690 
sec. b) t=61.710 sec. c) t=61.730 sec. d) t=61.750 sec. 

Figure 9. 

Surfaces of spatial correlation of fibrillation vs separation vector for exper- 
iment #3. Correlations were computed from 1.0 seconds of data at each 
sensor of a 22 x 23 electrode array, a) 1.0 seconds to 2.0 seconds after VF 
was induced, b) 2.0 - 3.0 seconds, c) 3.0 - 4.0 seconds, d) 60.0 - 61.0 seconds. 

Figure 10. 

Curves of spatial correlation vs magnitude of separation for experiment #3. 
The radial distance from the origin of Figures 8a-d) to each grid point is the 
magnitude of separation. See Table 2 for the correlation lengths correspond- 
ing to each curve. 

Figure 11. 

Surfaces of spatial correlation vs separation vector for experiment #4. a) 
1.0 - 2.0 seconds, b) 2.0 - 3.0 seconds, c) 3.0 - 4.0 seconds, d) 60.0 - 61.0 
seconds. 

Figure 12. 

Curves of spatial correlation vs magnitude of separation for experiment #4. 
See Table 2 for the corresponding correlation lengths. 

Figure 13. 

Mean (± standard deviation) correlation lengths over the five experiments 
at various stages of VF. The difference between correlation lengths at 60-61 
seconds and those at 3-4 seconds is statistically significant with p < 0.05. 
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